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We measure the cosmic microwave background (CMB) power spectrum on angular scales t ~ 
30 — 200 (1° — 6°) from the QMASK map, which combines the data from the QMAP and Saskatoon 
experiments. Since the accuracy of recent measurements leftward of the first acoustic peak is limited 
by sample- variance, the large area of the QMASK map (648 square degrees) allows us to place among 
the sharpest constraints to date in this range, in good agreement with BOOMERanG and (on the 
largest scales) COBE/DMR. By band-pass-filtering the QMAP and Saskatoon maps, we are able to 
spatially compare them scale-by-scale to check for beam- and pointing-related systematic errors. 



I. INTRODUCTION 

After the discovery of large-scale Cosmic Microwave 
Background (CMB) fluctuations by the COBE satellite 
, experimental groups have forged ahead to probe ever 
smaller scales. Now that TOCO ||,|, BOOMERanG |] 
and Maxima g have convincingly measured the location 
and height of the first acoustic peak, attention is shifting 
to still smaller scales to resolve outstanding theoretical 
questions. For instance, the structure of the second and 
third peaks constrains the cosmic matter budget^]. How- 
ever, it remains important to improve measurements on 
larger angular scales as well, both for measuring cosmo- 
logical parameters and to cross-validate different exper- 
iments against potential systematic errors. This is the 
goal of the present paper. 

Since the accuracy of recent measurements leftward 
of the first acoustic peak is limited by sample-variance 
rather than instrumental noise, we will use the largest 
area CMB map publicly available to date with degree 
scale angular resolution. This map, nicknamed QMASK 
H , is shown in Figure [l] and combines the data from the 
QMAP and Saskatoon ||-|5) (hereafter SASK 

sometimes) experiments into a 648 square degree map 
around the North Celestial Pole. This map has been 
extensively tested for systematic errors Q , with the con- 
clusion that the QMAP and Saskatoon experiments agree 
well overall, as well as analyzed for non-Gaussianity 
|^9|,^0|. However, for the present power spectrum anal- 
ysis, it is important to perform additional systematic 
tests to see if there is evidence of scale-dependent prob- 
lems in any of the maps. In particular, pointing prob- 
lems, beam uncertainties, sampling and pixelization ef- 
fects can smear the maps in a way that changes the shape 
of the power spectrum, suppressing small-scale fluctua- 
tions (and under some circumstances boosting the fluc- 
tuations as well). 




FIG. 1. Wiener-filtered QMASK map combining the QMAP 
and Saskatoon experiments. The CMB temperature is shown in 
coordinates where the north celestial pole is at the center of the 
dashed circle of 16° diameter, with R.A. being zero at the top and 
increasing clockwise. This map differs from the one published in 
Ji)| by the erasing of QMAP information for £ > 200 described in 
the text. 

The rest of this paper is organized as follows. In Sec- 
tion [n], we present a technique for erasing this type of 
unreliable information, and apply it to produce a new 
combined QMASK map where all the statistical weight 
on small scales comes from Saskatoon. We compute the 



power spectrum of this combined map in Section III 



* After this paper was first submitted, accurate new 
measurements of these higher peaks were reported by the 
Boomerang [0, DASI Q and Maxima || teams, so we have 
added comparisons with these results below. 



We discuss systematic errors in section IV, finding good 
agreement but a hint of suppressed QMAP power for 
I > 200, and conclude with a rather conservatively cut 
data set that we believe to be reliable as a starting point 
for cosmological parameter analysis. 
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II. COMBINING THE SASK AND QMAP 
EXPERIMENTS 



As will be discussed in more detail in Section [V 



point- 
ing problems, beam uncertainties, pixelization and sam- 
pling effects can suppress small-scale fluctuations in the 
CMB power spectrum. The QMAP results may there- 
fore only be valid on large angular scales |10| - [l2| . Both 
pointing inaccuracies and pixelization effects could have 
effectively smoothed the QMAP map, suppressing small- 
scale power. Flight 1 of QMAP [jlO| was sampled at a rel- 
atively low rate, causing the effective beam shapes to be 
elongated along the scan direction. Since this ellipticity 
was not modeled in the mapmaking algorithm jl2|, the 
resulting smoothing would again be expected to suppress 
small-scale power. Miller et al. |ll| present a detailed 
discussion of these issues, and conclude that the QMAP 
measurements are likely to be unreliable for £ > 200. 
Although attempts can be made to model and correct 
for some of these effects, an accurate treatment of the 
undersampling problem in particular would be way be- 
yond the scope of the present paper, requiring the entire 
maps to be regenerated from the time-ordered data with 
an order of magnitude more pixels to be able to resolve 
the beam ellipticity (the ellipticity is not uniform in di- 
rection, since the sky rotated relative to the scan axis 
during the flights). Such an analysis would hardly be 
worthwhile anyway, since the strength of QMAP com- 
pared to subsequent experiments lies on large scales, not 
on small scales. 

For these reasons, we adopt a more conservative ap- 
proach, combining the QMAP and Saskatoon maps in 
such a way that the small-scale QMAP information can 
be optionally erased as a precaution. By this we do not 
mean removing the signal (smoothing the map), which 
would just lead to further underestimation of the true 
power. Rather, we mean removing the information, i.e., 
doing something that causes subsequent analysis steps 
(like combining with Saskatoon or measuring the power 
spectrum) to give negligible statistical weight to the 
small-scale QMAP signal. We achieve this by creating a 
random map with very large small-scale noise and adding 
it to the QMAP map, modifying its noise covariance ma- 
trix N accordingly. 

In practice, we start by generating a white noise map 
x white which has the following properties: it covers the 
same sky region as QMAP, and each pixel temperature is 
drawn independently from a Gaussian distribution with 
zero mean and standard deviation a, giving it a noise 
covariance matrix H w hite = & 2 I. This makes its angu- 
lar power spectrum Ci independent of £. We then apply 
the Laplace operator V 2 to the mock map. Since it is 
pixelized on a square grid, we do this in practice by mul- 
tiplying by a matrix L that subtracts each pixel from the 
average of its four nearest neighbors. The transformed 
map Xbiuo = Lx w hitc thereby obtains a very blue power 
spectrum Ci since Laplace transformation corre- 



sponds to multiplying by £(£ + 1) in the Fourier (multi- 
pole) domain. We choose the normalization factor a such 
that the blue noise starts dominating the noise and sig- 
nal of the QMAP map around £ ~ 200. Since the CMB 
power falls off as Ci °S £~~ 2 , the result is that the added 
noise is negligible for I < 200 and dominates completely 
for I > 300. Finally, we add this blue noise map to the 
QMAP map, obtaining 



XQMAP + Xblu 



'QMAP 



->blu 



XQMAP + Lx w hite, (1) 
= SqmAP + cr 2 LL*. (2) 



To quantitatively assess the effect of this procedure, 
we perform a series of numerical experiments where we 
rescale the noise level by a factor p = 0, 1,3, 5, 7, 1000. 
This corresponds to multiplying Xbiuc by p and multiply- 
ing Sbiuc by p 2 . p — means adding no noise, i.e., that 
x M retains all QMAP information; p — oo (or p = 1000 
in practice) means that x new contains essentially pure 
noise. Below we will see that the choice p = 3 erases in- 
formation on small scales I > 200 very well, while doing 
little harm on larger scales. After erasing with p = 3, we 
combine the resulting maps with the Saskatoon data as 
in H . The result is shown in Figure [lj and looks almost 
unchanged (compare Figure 4 in 9]) since Wiener filter- 
ing suppresses noisy modes and the smallest scales were 
noisy to start with. 



III. THE ANGULAR POWER SPECTRUM 

In this section, we compute the angular power spec- 
trum of the QMASK map produced in the previous sec- 
tion, shown in Figure [j]. It contains 6495 pixels and 
covers a 648 square degree sky region. We calculate 
the angular power sp ectrum using the quadratic esti- 
mator method of |]l^ , ^o| , implemented as described in 
|§3>|§. This method involves the following steps: (i) S/N 
compression of data and relevant matrices by omitting 
Karhunen-Loeve (KL) eigenmodes with very low signal- 
to-noise ratio, (ii) computation of Fisher matrix and raw 
quadratic estimators, (iii) decorrelation of data points. 
We compute the power in 20 bands from I = 2 to 400 of 
width A£ = 20, which takes about a week on a worksta- 
tion. We repeated this calculation in its entirity for the 
following values of the p-parameter: (no information 
erased), 1, 3, 5, 7 and 1000 (Saskatoon only, QMAP fully 
erased), as shown in Figure 0. 

In the progression of power spectra shown in this 
figure, the ^-values beyond which QMAP is effectively 
erased is seen to shift to the left as p increases. For 
p = 3, this scale is seen to be of order 200 in the sense 
that the leftmost points are almost indistinguishable from 
the p = case whereas those with decent signal-to-noise 
for £ > 200 (say points number 10 and 11) are simi- 
lar to the p = 1000 case. Adopting this cutoff scale as 
our baseline calculation (we further discuss this choice 
in the next section), we then average these rather noisy 
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p — 3 measurements into six (still uncorrected) mea- 
surements listed in Table 2. The first four probe angular 
scales where the above-mentioned systematics are likely 
to be negligible, incorporating the first 9 £-bands (up to 
I = 180), and are shown in Figure [| The horizontal 
bars show the mean and rms width of the corresponding 
window functions. 
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FIG. 2. Angular power spectrum ST = [i(t + l)Ce/2n] 1 / 2 (un- 
correlated) in 20 bands of CMB anisotropy from the combined 
QMASK data in case of p = 0(no additional noise at all, just 
simple combination of SASK and QMAP),1, 3, 5, 7 and 1000 (no 
QMAP information, only SASK information). For comparison, we 
also plot the a recent "concordance" model [|l7j and the power 
measurements from MAXIMA and BOOMERanG. 
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bars are uncorrelated between the twenty measurements, but do 
not include an overall calibration uncertainty of 10% for ST. 
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FIG. 3. Uncorrelated measurements of the CMB power spec- 
trum ST = [£{l + 1)Q/2tt] 1 / 2 from the combined QMASK 
data. For comparison, we also plot the measurements from 
COBE/DMR, MAXIMA, DASI and BOOMERanG. These error 
bars do not include calibration uncertainties of 10% (QMASK), 
10% (BOOMERanG), 4% (DASI) and 4% (MAXIMA). As de- 
scribed in the text, we suggest using only the first four points 
(£ < 200), not the two dashed ones, for cosmological model con- 
straints. 
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STj [nK*] 


40 ± 16 
79 ± 13 
125 ± 19 
178 ± 15 


1308 ± 452 
1899 ±410 
2436 ± 471 
4295 ± 869 


259 ± 46 
335 ± 58 


4265 ± 726 
2596 ± 1795 



Table 1. The power spectra ST = {1(1 + l)C e /2n] 1 / 2 from the 
QMASK map in two cases:p = and p = 3. The tabulated error 



Table 2. The power spectrum ST = \t(t + l)C l /2ir} 1 / 2 from the 
QMASK map. The tabulated error bars are uncorrelated between 
the six measurements, but do not include an overall calibration 
uncertainty of 10% for ST. We recommend using only the first four 
for cosmological model constraints. 



IV. DISCUSSION 

We have measured the CMB power spectrum from the 
combined QMAP and Saskatoon data sets, obtaining sig- 
nificant detections in the range £ ~ 30 — 300. The key 
question that we need to address in this section is how 
reliable these measurements are. 

The "usual suspects" as far as systematics and non- 
CMB contamination are concerned tend to have power 
spectra that are either redder or bluer than that of the 
CMB, and therefore naturally split into two categories: 
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• Problems on large angular scales can be caused 
by contamination from diffuse foregrounds (syn- 
chrotron, free- free and spinning dust emission), and 
systematic errors related to atmospheric contami- 
nation, scan-synchronous offsets, atmospheric con- 
tamination, etc. 

• Problems on small angular scales can be caused 
by contamination from point source foregrounds 
and systematic errors related to pointing problems, 
beam uncertainties, pixclization and sampling ef- 
fects. 

The foreground contaminatio n pq j has been previously 
quantified for both Saskatoon ]27] and QMAP p8| , by 
cross-correlating the maps in question with a foreground 
templates tracing synchrotron, dust and free-free emis- 
sion, point sources, etc. The conclusion was that fore- 
grounds contribute at most a few percent to the angu- 
lar power spectrum reported here, mainly on the largest 
angular scales. The errors reported on the power spec- 
trum assume that the underlying CMB signal is Gaus- 
sian, which is supported by two recent Gaussianity anal- 
yses of the QMASK map ||& 

As opposed to foregrounds, which are true features on 
the microwave sky, the remaining problems listed above 
are experiment-specific effects which would be expected 
to affect QMAP and Saskatoon differently. This provides 
us with a powerful tool with which to test for their pres- 
ence, which we will now employ: comparing the QMASK 
and Saskatoon maps where they overlap. 

Previous comparisons between CMB data sets have 
been done either spatially or in terms of power spectra 
(discarding phase information). Since two inconsistent 
maps can have identical power spectra, one should be 
able to obtain still stronger tests for systematic errors by 
comparing power with phase information. For instance, 
one could imagine band-pass filtering the two maps to 
retain only a particular range of multipolcs I, and then 
testing whether these two filtered maps were consistent. 
We will now describe a simpler way of implementing a 
test in this spirit. 

A. A method for scale-by-scale comparison of two 
maps 

Numerous map comparisons have been performed with 
the "null-buster" test |l6) 

z'N-^SN^z-trfN^S} 

v = 1 m , 3 

[2 trjN-iSN-iS}] 172 

where v can be interpreted as the number of "sigmas" 
at which the difference map z is inconsistent with pure 
noise. If the two maps are stored in vectors xi and X2 
and have noise covariance matrices Nj and N 2 , then a 
weighted difference map z = Xx — rx 2 will have noise 



covariance N = Ni + r 2 N 2 . The matrix S tells the test 
which modes (linear combinations of the pixels z) to pay 
most attention to, and can be chosen arbitrarily. The 
choice S = N gives a standard x 2 -test. It can be shown 
]I6[ that the null hypothesis that z is pure noise (that 
(zz*) = N) is ruled out with maximal significance on 
average if S is chosen to be the covariance of the expected 
signal in the map, i.e., S = (zz f ) — N. In our case, we 
choose S to be the CMB covariance matrix corresponding 
to a power spectrum 

^rji J 1 /^K if £ G [£min:^max\: 

1 [0 /iK otherwise. ' ^ ' 

ensuring that the test only uses information in the mul- 
tipole interval [Imim&max]- The overall normalization of 
S is irrelevant, since it cancels out in equation (|J). 

B. Results of comparing QMAP and Saskatoon 
scale- by-scale 

We will now use the method described above to com- 
pare the QMAP and Saskatoon data scale-by-scale, i.e., 
in different multipole intervals. 

The QMASK map has been shown to be inconsistent 
with noise at the 62a level [§. In the region where the 
QMAP and Saskatoon maps overlap, they were found 
to detect signal at 40c and 21a, respectively, while the 
difference map was consistent with pure noise. Which 
angular scales are contributing most of this information, 
and how well do the two maps agree scale-by-scale? 




0.001 0.01 0.1 1 10 100 

Relative normalization r 



FIG. 4. Comparison of the QMAP and Saskatoon experi- 
ments on different angular scales, corresponding to the multipole 
ranges shown in square brackets. The curves show the number 
of standard deviations ("sigmas") at which the difference map 
xqmap — 7*xsask is inconsistent with mere noise. Note that this 
is only for the spatial region observed by both experiments. 

To answer these questions, Figure |] shows the result of 
comparing QMAP with Saskatoon in the four multipole 
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intervals [2,100], [101,200], [201,300], and [301,400]. 
QMAP is seen to detect signal in the overlap region at 
47(7, 18cr, 12(7 and 6(7, respectively, whereas the corre- 
sponding numbers for Saskatoon are 18(7, 14(7, 11(7 and 
7(7. In other words, QMAP dominates on large scales, 
whereas Saskatoon has at least comparable information 
content on small scales because of superior angular reso- 
lution. 

Although QMAP and Saskatoon both detect signifi- 
cant CMB signal in all four bands, this signal is seen to 
be common to both maps since the difference maps z for 
r = 1 are consistent with noise. Furthermore, there is 
no evidence of relative calibration errors for the [2, 100] 
or [101, 200] bands, since the minima of these two curves 
are at r « 1. However, the situation is less clear on 
smaller scales: the best-fit amplitude of QMAP is only 
63% of the amplitude of Saskatoon for £ G [201, 300], and 
even lower for £ G [301, 400]. None of these departures of 
the minimum from r = 1 are statistically significant — 
we cannot determine whether this is a a problem or not 
simply because the amount of information in the maps 
drops sharply on small scales where detector noise and 
beam dilution become important. However, whereas the 
Saskatoon information was extracted from highly over- 
sampled calculations of the relevant beam patterns on 
the sky, and should be reliable on small scales, there are 
a number of reasons why the QMAP data may only be 
valid on larger [i < 200) angular scales flC)|-|l2]| : 

1. The QMAP pointing solution was only accurate to 
this level El|, and small residual pointing uncer- 
tainties could have effectively smoothed the map, 
suppressing power for £ > 200 relative to Saska- 
toon just as Figure ^ indicates. 

2. The QMAP maps were generated by subdividing 
the sky into square pixels of side 9 = 0.3125°, and 
the effect of this pixelization may well become im- 
portant on angular scales substantially exceeding 
£ ~ 1/6 k, 200, suppressing power on these scales. 

3. Flight 1 of QMAP @, which dominates the sky 
coverage in Figure [l] was sampled at a relatively low 
rate, causing the effective beam shapes to be elon- 
gated along the scan direction. Since this ellipticity 
was not modeled in the mapmaking algorithm |l2) , 
the resulting smoothing would again be expected 
to suppress power on scales £ > 200. 

Miller et al. |l8| present a detailed discussion of these 
issues, and conclude that the QMAP measurements are 
likely to be unreliable for £ > 200. This justifies our p = 
3 choice in Section 2, effectively discarding this suspect 
small-scale information from QMAP. 

C. Conclusions 

We have measured the CMB power spectrum from 
the combined QMAP and Saskatoon data sets, and per- 



formed a number of tests for potential systematic errors. 
On large angular-scales, £ < 200, our results appear very 
solid: foreground contamination has been quantified to 
contribute no more than a few percent to the power spec- 
trum, and we find beautiful internal consistency between 
the QMAP and Saskatoon components of our map, both 
in terms of power amplitude and in terms of spatial phase 
information. On smaller scales I > 200, there are good a 
priori reasons to discard the QMAP information, and we 
have therefore done so. Our internal scale-by-scale com- 
parison between QMAP and Saskatoon independently 
suggests that there may be systematic problems on these 
scales. Since this unfortunately leaves us with no inde- 
pendent way of validating our Saskatoon results in this 
small-scale regime, e.g., testing our Saskatoon deconvolu- 
tion procedure, we strongly recommend the conservative 
approach of using only the first four band power mea- 
surements we have reported in table 2 and Figure 3 - 
this is also where QMASK is most sensitive relative to 
other experiments. These four measurements, spanning 
the range £ ~ 30 — 200, provide among the sharpest and 
best tested constraints to date on these scales, provide a 
good starting point for constraining cosmological models. 

D. Comparison with other experiments 

Figure || shows that our results agree well with the 
"concordance" model of flffl. They are also consis- 
tent with BOOMERanG ||7 DASI [Q and Maxima j| 
once calibration uncertainties are taken into account. 
The QMAP and Saskatoon calibration uncertainties are 
6% — 10% and 10 %, resp ectively (we have corrected the 
original ST results |10|,|l4| by a factor 1.05 using the latest 
Cassiopeia A data [|23[~as in (24)). Although the uncor- 
rected components of this may average down somewhat 
when the two maps are combined, we quote a 10% un- 
certainty on our result to be conservative. Our results 
are also consistent with those obtained from QMAP and 
Saskatoon alone. Our leftmost data point agrees well 
with the last point measured from COBE/DMR ||] by 
II- 

The original BOOMERanG98 power spectrum was 
based on a map area of 436 square degrees jl, so one 
might expect our error bars on large scales to be a factor 
(648/436) 1/2 « 1.2 smaller. Our actual error bars on STf 
are only about 10% smaller on the best QMASK scales af- 
ter adjusting for bandwidth differences, which is because 
of BOOMERanG's lower noise levels (the scan strat- 
egy and 1 / /-noise of QMAP introduced a non- negligible 
amount on noise even on the largest scales). We note 
that the substantial reduction in error bars relative to 
the original Saskatoon analysis |l4|] is due to additional 
information not only from QMAP, but from Saskatoon 
as well. This is because our present method extracts all 
the information present, whereas that employed in |l4] ] 
was limited to information along radial scans, not using 
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phase information between scans. 

In conclusion, we have measured the CMB power spec- 
trum on angular scales t ~ 30 — 200 from the QMASK 
map, placing among sharpest and best tested constraints 
to date on the shape of the CMB power spectrum as it 
rises towards the first acoustic peak. Our window func- 
tions, the combined map and its noise covariance matrix 
are available at 

www.hep.upenn.edu/ ~ xuyz / qmask.html. 
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